arXiv:astro-ph/0503050vl 2 Mar 2005 


X-ray Outburst in Mira A 


Margarita Karovska 1 , Eric Schlegel 1 , Warren Hack 2 , John Raymond 1 , and Brian E. Wood 3 * 
Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 


ABSTRACT 

We report here the Chandra ACIS-S detection of a bright soft X-ray transient in the 
Mira AB interacting symbiotic-like binary. We resolved the system for the first time 
in the X-rays. Using Chandra and HST images we determined that the unprecedented 
outburst is likely associated with the cool AGB star (Mira A), the prototype of Mira- 
type variables. X-rays have never before been detected from an AGB star, and the 
recent activity signals that the system is undergoing dramatic changes. The total X-ray 
luminosity of the system is several times higher than the luminosity estimated using 
previous XMM and ROSAT observations. The outburst may be caused by a giant flare 
in Mira A associated with a mass ejection or a jet, and may have long term consequences 


on the system. 

Subject headings: binaries: symbiotic— stars: 

Individual (Mira AB) — stars: winds, outflows 

1. Introduction 

Mira AB is an interacting binary system 
composed of an aging cool giant (Mira A) los¬ 
ing mass at a rate of ~1CU 7 M so i/yr (Bowers 
& Knapp, 1988), and an accreting compan¬ 
ion (Mira B) about ~70 AU (~0.6”) away. 
This detached binary is one of very few wind 
accreting systems that has been spatially re¬ 
solved and for which the energy distribution 
of both components can be determined unam¬ 
biguously (Karovska et al, 1997). Studies of 
Mira AB wind accretion and mass transfer at 
wavelengths ranging from X-rays to radio pro- 
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activity — stars: AGB and post-AGB — stars: 
X-Rays: general 

vide a basis for understanding wind accretion 
processes in many other astronomical systems 
that currently cannot be resolved. There¬ 
fore this nearest symbiotic-like system offers a 
test-bed for detailed studies of wind accretion 
processes and accretion theory (Livio 1988). 

In 1995, using the HST FOC, we resolved 
the system for the first time spatially and 
spectrally at UV and optical wavelengths and 
studied the interacting components individ¬ 
ually (Karovska et al. 1997). In the past 
few years we have been witnessing changes in 
the spectral energy distribution (SED) of Mira 
AB, especially in UV (Karovska et al. 1997; 
Wood, Karovska, &Hack 2001). In addition to 
the general fading of the accretion luminosity, 
another baffling development was the appear¬ 
ance of a forest of Lya-fluoresced H 2 emission 
lines, which dominated the HST spectra in 
1999, despite not being seen at all in the 1995 
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spectra or by IUE (Wood, Karovska, & Hack 
2001; Wood, Karovska, & Raymond, 2002). 
A similar drop in the accretion luminosity and 
appearance of a set of Lya-fluoresced H 2 emis¬ 
sion lines were also seen in the FUSE spectra 
in 2001 (Wood & Karovska 2004). Variable 
mass loss from Mira A could have caused the 
changes in the UV flux of Mira B observed in 
the past 10 years signaling that the system is 
undergoing dramatic changes. 

In 1993, the ROSAT observation of Mira 
AB resulted in the first unambiguous de¬ 
tection of X-ray emission from Mira AB 
(Karovska at al. 1996). This observation 
resolved the contradicting results from the 
analysis of the EINSTEIN observation; Jura 
& Helfand (1984) marginally detected an X- 
ray source, while Maggio et al. (1990) set an 
upper limit of f x <lA xl0~ 13 erg s -1 cm -2 . 
The ROSAT X-ray luminosity of the Mira AB 
system was estimated ~10 29 erg s^ 1 (Karovska 
et al. 1996) which is similar to the luminosity 
estimated from the XMM observations car¬ 
ried out about ten years later (Kastner & 
Soker 2004b). 

In this Letter we report the initial results 
from our recent Chandra observations and 
from coordinated HST and ground-based ob¬ 
servations of Mira AB, including our discovery 
of a soft X-ray outburst in Mira A. 

2. Observations and Analysis 

On 2003 December 6 we carried out a 70 
Ksec pointed Chandra observation of Mira AB 
using the ACIS-S instrument (Weisskopf et al. 
2002). The source was placed at the nominal 
aim point of ACIS-S, on CCD S3. We ana¬ 
lyzed the Chandra observations using Cl AO 
data reduction and analysis routines. 1 The 
source spectrum was extracted using an aper- 

1 CIAO is the Chandra Interactive Analysis of 
Observation’s data analyses system package 
(http://cxc.harvard.edu/ciao) 


ture of 12 7/ radius centered on Mira AB; the 
background came from a surrounding annulus 
of inner radius 14 7/ and outer radius 30A Fig¬ 
ure 1 shows the Chandra spectrum and the 
best spectral fit (described in detail in §3.2). 
The source contained 5115 counts for a count 
rate of 0.073T0.001 cts s -1 . The background 
rate was 0.0007 cts s -1 . In the 0.2-0.7 keV 
band, the count rate was 0.066T0.001 cts s -1 
(background: 0.0003 cts s _1 ); in the 0.7-4 
keV band, the detected rate was 0.008 cts s -1 
(background: 0.0004 cts s -1 ). The combina¬ 
tion of the declining effective area below 0.5 
keV plus the excess absorption layer (Chan¬ 
dra Proposer’s Guide 2 ) ensure minimal or zero 
photon pile-up. Pushing all of the model pa¬ 
rameters to their most pile-up-favorable ex¬ 
tremes yields a pileup estimate of less than 
4%. 

We detected several thousand counts below 
1 keV associated with a new bright soft source 
in the system, which was not reported a few 
months before by XMM (Kastner & Soker, 
2004b), or by ROSAT in 1993 (Karovska 
et al., 1996). The high-energy component 
(>lkeV) is similar in appearance to the qui¬ 
escent XMM and ROSAT spectra. However, 
a detailed comparison shows that a clear evo¬ 
lution has occurred in this portion of the spec¬ 
trum (~1.2-1.8 keV) as well. The count rates 
(see below), if folded through the XMM or 
ROSAT effective areas, would have been eas¬ 
ily detected in either instrument. Further¬ 
more, the spectrum is dominated by the soft 
component in the 0.2-0.7 keV band and the 
spectral behavior of the ROSAT and XMM 
spectra is significantly harder. 

Following the Chandra X-ray outburst de¬ 
tection, we requested and obtained Direc¬ 
tor’s time on Chandra and HST. Chandra 
HRC/LETG observations (40 Ksec) were car¬ 
ried about two months after the initial Chan- 

2 http://asc. harvard.edu/proposer/POG/iridex.html 
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dra observations (on January 11 2004). We 
detected signal only in the zeroth order of the 
spectrum which showed that the brightness 
of the soft source had dropped significantly. 
This observation detected ~80 counts. Were 
the soft ACIS spectral component (§3.2) still 
present with the same absorbed flux at the 
time of the LETG observation, >250 counts 
would have been detected in the zero-order 
image. However, the detected counts could 
be due to the hard source as well, or to both 
sources. Further modeling using the HRC and 
the ACIS observations is in progress. 

We carried out follow-up HST STIS CCD 
and MAMA observations about two months 
after the initial Chandra observations (2004 
February 2 and February 16). We resolved 
the system in the near-UV and obtained in¬ 
dividual spectra of both components. Mul- 
tiDrizzle techniques (Koekemoer at al. 2002) 
were used to combine multiple CCD images 
recorded using the F28x50OII filter centered 
on the [Oil] 3729 A line. The images resolved 
the system and detected extended emission. 
The STIS spectra show a significant increase 
in the Mg h & k lines emission in Mira A 
and in Mira B when compared to the 1999 
STIS observations. Velocity shifts of the or¬ 
der of 100 km/s were also detected in Mira 
A’s Mg h & k lines. A similar increase of line 
emission was detected in the optical spectra 
(H q and OIIIA5007A lines) obtained at the 
Oak Ridge observatory. In addition, broad¬ 
ening and complex structure were observed 
in the H a during the month following the X- 
ray outburst, indicating velocities of ~ 100 
km/s. Further details of the analysis and 
the results of the Chandra, HST, and ground- 
based observations will be discussed elsewhere 
(Karovska et al ., 2005, in prep.). 


3. Results 

3.1. X-ray and near-UV Imaging 

The ACIS-S raw image of Mira AB binned 
at the detector 0.495” pixel resolution showed 
an elongation along the Mira AB system axis 
(Figure 2a). We made soft (< 0.7 keV) and 
hard (0.7-2 keV) images of Mira AB, filtered 
based on the major components detected in 
the ACIS-S spectrum. Figure 2b shows the 
image of the hard source (to the East) with an 
overlay of the contours of the soft component 
(to the West). The centroids of these sources 
are displaced by 

We attempted to do further spatial analysis 
using the HRC zeroth order image. The image 
appears extended and it is shifted ~0.3 arc- 
seconds from the soft ACIS image, in the di¬ 
rection toward the hard ACIS image. It could 
therefore be associated with Mira A or Mira 
B (or both). However, the low counts and 
lack of spectral resolution make it difficult to 
determine if there is a component associated 
with a remnant of the outburst. Although 
the HRC pixels (~0.14”) are indeed smaller 
than the ACIS pixels, the actual resolution of 
the HRC is in fact comparable to ACIS. Fur¬ 
thermore, there are additional artifacts in the 
HRC images that makes it very difficult decid¬ 
ing if the extension in the image to the South 
and an additional source 0.5” to the West are 
real. 

We explored the spatial extent of the emis¬ 
sion in the ACIS-S (0.3-4 keV) image of the 
system at 0.2” resolution using a new mul¬ 
tiscale deconvolution technique EMC2 (Esch 
et al. 2004) and Chandra PSF model. This 
was possible because the Chandra data in¬ 
clude information about the photon energies 
and positions which was used to obtain fil¬ 
tered images and carry out sub-pixel resolu¬ 
tion analysis. The Chandra PSF varies as a 
function of energy and off-axis angle; we car¬ 
ried out PSF simulations using CIAO software 
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and threads, including the ChaRT PSF sim¬ 
ulator 3 . The PSF simulation was carried out 
using information on the spectral distribution 
and off-axis location of the system as inputs to 
ChaRT. The EMC2 deconvolution technique, 
described in detail in Esch et al. 2004, was 
developed for low-count statistics data, and 
it provides error estimates in addition to the 
reconstructed images. 

The deconvolved image (color image in Fig¬ 
ure 3) shows two sources separated by ~0.6”. 
The pixel size is 0.1” (0.2 ACIS-S pixel size). 
This is the first image of an interacting bi¬ 
nary that has been spatially resolved at X- 
ray wavelengths. We note that the location 
of the brighter source in the deconvolved im¬ 
age (to the West) corresponds to the position 
of the soft source (as determined using filtered 
ACIS-S images), and the fainter source (to the 
East) corresponds to the hard-band ACIS-S 
image. 

We compared the X-ray images with the 
HST images of the Mira AB components ob¬ 
tained two months later, showing two sources 
separated by ~0.6”. In the HST images Mira 
A is to the west of Mira B. Figure 3 displays 
the overlay of the contours of the HST 3729 
A image on the Chandra deconvolved image, 
showing that the soft X-ray source is in the 
vicinity of the of Mira A. The soft source 
is therefore likely associated with the AGB 
star rather than with the accreting companion 
Mira B. When the X-ray image is compared to 
the HST image of Mira A and Mira B, we see 
slightly larger separation of the X-ray compo¬ 
nents (~0.1”) and possible small rotation of 
the axis between the sources. This difference 
could be real, but also it could be due to un¬ 
certainties in the deconvolution, or to the lim¬ 
ited resolution in the Chandra data. It could 
also be a result of HST absolute positioning 
uncertainties. 

3 http://asc. harvard.edu/chart/index.html) 


The X-ray images of Mira A and Mira B ap¬ 
pear extended. Furthermore, the Chandra im¬ 
age shows a faint “bridge-like” feature extend¬ 
ing between the components. Similar “bridge¬ 
like” structure was detected in the 1995 HST 
images indicating possible mass flow between 
the components (Karovska et al. 1997). Al¬ 
though the features in the X-ray image are un¬ 
certain because of the low count level and/or 
possible PSF and deconvolution artifacts, we 
note that very similar structures can be seen 
in the overlayed HST contours. 

We determined that some of the extended 
structure in the HST image to the right of 
Mira A, in the opposite direction from Mira 
B, is due to an undocumented red leak of the 
F28x50OII filter. Following the HST obser¬ 
vations we discovered and confirmed the red 
leak in other very red sources. For example, 
F28x50OII filter images of CH Cyg, which 
contain a very red component, show similar 
structure, and we have detected similar struc¬ 
tured PSF in the [OIII] 5007 filter during the 
1999 observation of Mira AB. 

The evidence of extended structures in 
both components in both X-ray and UV im¬ 
ages is uncertain and require further analysis 
and modeling. In order to determine if there 
was a mass ejection during the outburst, it 
will be necessary to carry out further high- 
angular resolution monitoring including at X- 
ray and UV wavelengths. 

3.2. X-ray Spectroscopy 

The spectrum was binned to a minimum of 
20 counts per channel. The effective area was 
corrected for the energy- and time-dependent 
instrumental absorption and Xspec vll.O was 
used to fit the spectra (Arnaud 1996) as de¬ 
tailed below. The X-ray spectrum shows 
strong soft emission and potentially multi¬ 
ple harder components. Standard continuum 
models (bremsstrahlung, blackbody, power 
law) did not provide good fits to the soft com- 
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ponent. We describe the spectral fit below 
with 90% errors on the model parameters. 

We fit the low-energy portion of the spec¬ 
trum using an absorbed model spectrum con¬ 
sisting of a sum of gaussians, each with a 
fixed center and zero width to represent an 
un-resolved line. We included several C, N, 
and O lines. The line centers (in keV) are at 
C V 0.299, 0.304, 0.308; C VI 0.367; N VI 
0.319, 0.426, 0.431; C VI 0.435, 0.459; N VII 
0.500; O VII 0.561, 0.569, 0.574; and O VIII 
0.653). For each line, only the normalization 
was fit. The initial fit showed that the O line 
normalizations were all consistent with zero 
as were approximately half of the C and N 
normalizations; a subsequent fit used a model 
containing only the C and N lines. Figure 
1 shows the result with lines at C VI 0.367 
keV, N VI 0.426 keV, C VI 0.459 keV, and 
N VII 0.500 keV (line normalizations = 4.9, 
6.2, 50.8, and 33.1, respectively, in units of 
10~ 5 photons cm -2 s -1 ). Because of the lim¬ 
ited spectral resolution we cannot determine 
if there are any specific lines detected, or if 
the fit is unique. When we added a contin¬ 
uum component, its model normalization was 
consistent with zero. The apparently large 
residuals near 0.3 keV may be explained by 
the known problems with the corrections for 
the contaminant near the carbon edge (Mar¬ 
shall 2003 4 ). The best-fit spectrum is there¬ 
fore most easily explained as blended emission 
of C + N lines. 

The hard component was best-fit with a 
bremsstrahlung continuum of kT ~0.78^o23 
keV to mimic an optically thin thermal 
plasma, plus two zero-width gaussians rep¬ 
resenting emission lines at 1.02^{j'[j(j and 
1.36 ^q'q 4 keV, all absorbed by a column 
Nh ~9.0^q 42x 10 22 cnb 2 . Other models 
provided very poor fits including models 
containing additional gaussians at the posi¬ 


tions of expected line emission. The fitted 
bremsstrahlung temperature is nearly identi¬ 
cal with the value determined from the XMM 
spectrum (Kastner & Soker 2004b). The fit¬ 
ted column is about a factor of two higher 
than the XMM-determined value. A specific 
error range is not shown for the XMM value, 
but the authors state “formal uncertainties 
are ^20%”. A 20% error range on the XMM 
value places it outside the 90% range of the 
fitted value determined here suggesting the 
column may have increased by a factor of 
~2 between the two observations. The hard 
gaussians have equivalent widths of 236^4 4 q 
and 124lgQ 4 eV, respectively, and we attribute 
them to Ne and Mg emission lines. We note, 
Ne IX and Ne X lines were also detected by 
the XMM (Kastner & Soker 2004), but the 
Mg line is not pronounced in the XMM spec¬ 
trum. These results are similar to the en¬ 
hanced Ka H-like Ne and He-like Mg lines 
observed in the ASCA spectrum of the sym¬ 
biotic CH Cyg showing jet activity (Ezuka et 
al. 2001), corresponding to log(T)=6.8. Our 
signal-to-noise is insufficient to fit for the He- 
like Si line (found in the CH Cyg spectrum). 
These high temperature lines may indicate 
shock heated emission or could be associated 
with the flare or a jet-like ejecta. 

We calculated the total flux from the 
system 5.6xlO -13 (0.2-4 keV), 5.5xlO” 13 
(soft component only in 0.2-0.7 keV), and 
2.7xl0~ 14 (hard component only in 0.7-4 
keV) erg s _1 cnr 2 for the bands defined. 
The unabsorbed fluxes are very sensitive to 
the adopted column. The fitted column to¬ 
ward the soft source is ~8.0xl0 19 cm 2 . This 
value, if adopted as the column toward the 
system, is not consistent with the column ob¬ 
tained from the ROSAT spectrum (Karovska 
et al. 1996) nor the XMM observation (Kast¬ 
ner & Soker 2004b). However, the fitted col¬ 
umn is really only a lower limit. The 90% 
contour is closer to ~2.8xl0 20 
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cm -2 . Within this upper limit, the column 
remains inconsistent with the XMM value, 
but is consistent with the column determined 
from the UV H 2 line spectroscopy (Wood et 
al. 2002). If we adopt the formal-fit col¬ 
umn from the soft spectrum, then the un¬ 
absorbed fluxes are 6.4xl0 -13 (0.2-4 keV), 
6.2x 10 -13 (0.2-0.7 keV), and 2.7xl0~ 14 (0.7- 
4 keV) erg s~ 4 cnh 2 . If we adopt the upper 
limit value for the column, the unabsorbed 
fluxes are 8.3xlO -13 (0.2-4 keV), 8.0xlO -13 
(0.2-0.7 keV), and 2.7xl0 -14 (0.7-4 keV) erg 
s^ 1 cm~ 2 . Adopting the fitted value of the 
column for the hard component, the unab¬ 
sorbed hard band flux is 9.9xl0~ 14 erg s -2 
cm -2 (0.7-4 keV). 

We emphasize that the appropriate Nh 
value could be completely different for the 
hard and soft sources since the emissions orig¬ 
inate from different places. It is conceivable 
that a persistent low energy X-ray source ex¬ 
ists, and that it was made visible by a tempo¬ 
rary drop in the neutral column density sur¬ 
rounding Mira A. It seems more likely, how¬ 
ever, that an X-ray flare ionized its surround¬ 
ings, producing the formal value of smaller 
column as a side effect. Taking note of the 
different columns toward the sources, the un¬ 
absorbed luminosity then falls into the range 
of 1.5-2 x 10 30 erg s -1 , for an assumed distance 
of 130 pc, and assuming the luminosity is the 
sum of the soft and hard components sepa¬ 
rately. The range results from the adopted 
value of the column toward the soft compo¬ 
nent. The total luminosity is about a factor 
of 2-4 higher than that observed by XMM and 
ROSAT. A detailed analysis and modeling of 
the Chandra data combined with the multi¬ 
wavelength observations of the outburst will 
be discussed elsewhere (Karovska et al, 2005, 
in prep.). 


4. Possible Causes and Consequences 
of the X-ray Outburst 

Detecting X-ray emission from Mira A is a 
very important result because X-ray emission 
has not been detected so far in single AGB 
stars, including Miras (e.g. Kastner & Soker 
2004a). The soft X-ray outburst in Mira A 
could be caused by a magnetic flare followed 
by a large mass ejection, analogous to the so¬ 
lar flares and Coronal Mass Ejections(CMEs). 
Soker and Kastner (2003) suggested that mag¬ 
netic reconnection events near the stellar sur¬ 
face should lead to localized, long-duration 
flares. The crude time scale for a flare might 
be estimated by multiplying the 1-2 day time 
scale for a large flare on an RS CVn star by 
the ratio of Mira A’s radius (few hundred solar 
radii) to that of the RS CVn star (4-15 solar 
radii), suggesting time scales from weeks to 
months, which is in accordance with the time 
scale of the X-ray event observed by Chandra. 

In the case of mass ejection we could ex¬ 
pect changes in the SED of both components 
on a time scale of months to years. As the 
ejected material cools down we expect a sig¬ 
nificant increase in dust formation in the sys¬ 
tem in 2004/2005. This is also the time scale 
on which we would expect a response of the 
accretion disk around Mira B to any distur¬ 
bance caused by an outburst on Mira A in 
late 2003, if the flow is propagating toward 
Mira B with speed of ~ hundred krn/s (as 
suggested by our preliminary analysis of the 
HST and ground-based spectroscopy). An in¬ 
creased outflow from Mira A resulting from 
the recent outburst could help renew the ac¬ 
cretion disk around the companion, following 
the low-state observed by HST and FUSE in 
1999-2001. 

The increased accretion rate into Mira B 
could also cause instabilities in the accretion 
disk and jet-like activity (e.g. Soker 2002) 
similar to that detected in several nearby 
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unresolved symbiotic systems (e.g. in CH 
Cyg, Corradi et al. 2001, R Aqr, Kellogg 
et al. 2001). An outburst in Mira B as¬ 
sociated with the accretion disk instabilities 
could be analogous to a dwarf nova outburst 
except that the much larger size of the wind- 
fed disk in this well-separated binary would 
increase the duration from about a week to 
many months. Further multiwavelength ob¬ 
servations are necessary to determine the na¬ 
ture and characteristics of instabilities in the 
system and/or its individual components that 
may have caused the outburst, and to under¬ 
stand the long-term consequences of this out¬ 
burst on the system. 
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Fig. 1.— Chandra ACIS-S spectrum of Mira 
AB (plus signs) fit with a combination of gaus- 
sians for the soft spectral component and a 
bremsstrahlung plus gaussians for the hard 
component (see §3.2). Residuals fall mostly 
near the C edge at 0.3 keV where the response 
matrix is known to have errors. 






Fig. 2.— Chandra images of Mira AB: (a) 
ACIS-S raw image of Mira AB filtered from 
0.3 to 2 keV. Mira A is toward the West (see 
Fig. 3); (b) contours of the ACIS-S soft im¬ 
age (0.3-0.7keV) (toward the West) overlayed 
on the image of the hard image (0.7-2 keV) 
(toward the East). 
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Fig. 3.— Chandra image of Mira B (left) and 
Mira A (right), separated by ~ 0.6”, with 
overlayed contours of the HST 3729 A image 
of the system. North is up, East is to the left. 
The apparent extended point-like structures 
in the HST contours of Mira A in the NW di¬ 
rection are PSF structures due to a red leak 
in the filter. 
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